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ABSTRACT 


In  this  work,  the  exterior- interior  coupling  problem  of  a cavity 
backed  aperture  in  a perfectly  conducting  infinite  sheet  is  considered. 
The  cavity  is  assumed  to  be  rectangular  and  contains  a perfectly  conduct 
ing  obstacle.  For  numerical  considerations  the  obstacle  is  taken  to  be 
a straight,  thin  wire,  oriented  perpendicular  to  one  of  the  cavity  walls 
The  problem  is  formulated  in  the  frequency  domain  with  an  elwT  time 
dependence. 

The  dyadic  formulation  of  this  vector  electromagnetic  boundary 
value  problem  is  given.  The  controversy  over  the  longitudinal  wave 
functions  and  their  contribution  to  field  dyads  is  resolved.  Specific 
attention  is  given  to  the  singularities  and  completeness  of  the  Green's 
dyads'  eigenfunction  expansions. 

Numerically  tractable  integral  equations  for  the  aperture 
electric  fields  are  obtained.  A summation  method  is  developed  for 
evaluation  of  otherwise  slowly  converging  eigenfunction  expansions  of 
potential  dyads  when  source  and  observation  points  become  close. 
Suggestions  for  future  extensions  of  this  work  are  discussed. 


CHAPTER  1 


r 


i. 


INTRODUCTION 

The  problem  considered  (Figure  1)  is  the  exterior-interior 
coupling  of  a cavity  backed  aperture  in  a perfectly  conducting  infinite 
sheet.  This  paper  extends  the  work  of  Seidel  [1977,  1978]  eliminating 
the  need  for  his  dipole  approximations  for  the  aperture  fields.  The 
cavity  considered  (Figure  1)  is  rectangular  with  "a"  less  than  or  equal 
to  "b",  has  perfectly  conducting  cavity  walls,  and  contains  a perfectly 
conducting  wire  obstacle.  In  the  general  formulation  (Chapter  4)  the 
obstacle  may  be  of  general  shape.  For  numerical  considerations  however 
it  is  taken  to  be  a thin  wire  parallel  to  two  of  the  cavity  walls.  The 
problem  is  formulated  in  the  frequency  domain  with  an  elti)T  time  depend- 
ence. In  order  to  obtain  the  summation  formulas  of  Chapter  6,  it  is 
assumed  that  the  frequency  of  operation  is  such  that  the  wave  number 
k times  the  cavity  dimension  "a"  is  less  than  it. 

Chapter  2 deals  with  the  dyadic  formulation  of  three  dimensional 
electromagnetic  boundary  value  problems.  In  particular  eigenfunction 
expansions  are  given  for  rectangular  cavity  dyads.  Completeness  of  the 
cavity  expansions  is  considered  as  well  as  their  singularities.  Finally 
the  contribution  of  the  longitudinal  vector  wave  functions  for  free 
space,  cavity,  and  waveguide  dyadic  expansions  is  examined. 

In  Chapter  3,  attention  is  given  to  preliminary  two  dimensional 
slot  diffraction  problems.  Sample  numerical  results,  illustrating  the 
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Figure  1.  Geometry  of  Problem 
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expected  field  behavior  for  the  three  dimensional  problem  (Figure  1)  are 
given.  For  one  incident  field  polarization  the  field's  edge  behavior 
has  been  built  into  the  numerical  solution.  This  illustrates  the  use 
of  the  edge  behavior  to  reduce  the  matrix  size  needed  to  solve  the 
discretized  problem.  The  incorporation  of  field  edge  behavior  in  the 
three  dimensional  problem  still  needs  to  be  studied. 

The  remaining  chapters  deal  with  the  numerical  formulation  of 
the  three  dimensional  coupling  problem  of  Figure  1.  In  Chapter  4 
techniques  of  Wilton  and  Dunaway  [1975]  are  extended  to  obtain  numeri- 
cally tractable  integral  equations  for  aperture  fields.  In  Chapter  5, 
the  method  of  moments  [Harrington  1968]  is  used  to  put  the  problem  in 
matrix  form.  In  the  analysis,  the  wire  obstacle  has  been  removed  since 
it  has  been  treated  by  Seidel  [1977,  1978].  Chapter  6 deals  with 
summing  otherwise  slowly  convergent  eigenfunction  expansions  of  poten- 
tial dyads  needed  for  carrying  out  the  numerical  procedure  of  Chapter  5. 
Finally,  Chapter  7 gives  conclusions  as  well  as  several  suggestions  for 
extension  of  this  work. 

Background 

In  recent  years  problems  for  which  a wire  couples  to  an  aperture 
in  an  infinite  perfectly  conducting  sheet  have  been  considered  [Butler, 
and  Umashankar  1976;  Seidel  and  Butler  1976].  However  the  image  theory 
needed  to  apply  the  techniques  put  forth  in  the  above  papers  to  the 
cavity  backed  aperture  problem  with  obstacle  yields  a three  dimensional 
array  of  images  [Seidel  1977,  p.  22].  Butler,  Rahmat-Samii  and 
Mittra  [1978]  provide  a good  review  of  aperture  literature. 
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An  alternate  method  for  formulation  of  the  problem  utilizes 
dyadic  Green's  functions.  An  expansion  for  the  electric  field  cavity 
dyad  was  found  by  Weyl  [1913,  1915],  Teichmann  [1952]  found  that  this 
expansion  was  incomplete  and  provided  a correction  [Teichmann  and 
Wigner  1953].  The  dyad  for  the  magnetic  vector  potential  for  the 
rectangular  cavity  is  formulated  in  Morse  and  Feshback  [1953,  pp.  1849- 
1851].  More  recently  these  dyads  have  been  formulated  by  Tai  and 


r 


Rozenfeld  [1976]  and  Rahmat-Samii  [1975].  However  there  has  been  some 
controversy  over  the  effect  of  the  longitudinal  £ vector  wave  function 
contribution  to  these  field  dyads.  This  point  is  covered  in  detail  in 
Chapter  2 and  the  difficulties  resolved.  (A  review  of  the  appropriate 
literature  is  given  in  Chapter  2.) 

Seidel  [1977,  1978]  has  solved  the  problem  (Figure  1)  under  the 
small  aperture  approximation  by  the  use  of  dipole  moments.  The  major 
contributions  of  the  present  work  are  the  removal  of  these  dipole 
approximations,  resolution  of  the  controversy  over  the  longitudinal 
wave  functions,  the  derivation  of  numerically  tractable  integral 
equations  for  the  aperture  fields,  and  the  summation  of  the  potential 
dyads'  expansions  (Chapter  6]  when  the  distance  between  source  and 
observation  points  is  small. 


1 1 
* 1 


r 

r> 


CHAPTER  2 


VECTOR  FIELDS  AND  DYADIC  FORMULATION 

The  problem  considered  is  the  computation  of  steady  state  (eltoT) 
electromagnetic  fields  in  a region  with  known  sources.  The  region's 
boundary  (Figure  1)  is  assumed  to  be  one  of  the  following:  the  sphere 
at  infinity,  or  a perfectly  conducting  surface  Sc  with  apertures,  or 
part  of  the  sphere  at  infinity  plus  a perfectly  conducting  surface  Sc 
with  apertures. 

Maxwell's  equations  for  steady  state  electric  sources  in  a 
homogeneous,  linear,  isotropic  medium  may  be  written  as: 


V x E = H 

(2.1a) 

V x H = k2  H + j 

(2.1b) 

~ E 

E = -4 — 

(2.1c) 

-ltoy 

* P 

V • E = “r“ — 

(2. Id) 

-icjgu 

V • H = 0 

(2. le) 

The  magnetic  source  analogues  are 

V x H = E (2.2a) 

VxE=k2H+M  (2.2b) 
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The  above  equations  yield  the  vector  wave  equations 


V x V x E - k2  E = J (2. 

VxVxH-k2H=VxJ  (2. 

and 

VxVxE-k2E=VxM  (2. 

7 x V x H - k2  H = M (2. 


for  electric  and  magnetic  source  excitations,  respectively. 

A vector  analogue  of  Green's  second  identity  [Stratton  1941, 
p.  250]  will  be  used  to  invert  these  wave  equations,  viz: 


VxVxP  - P • V x V x Q)dV 


(PxVxQ-QxVxP) 


n dS  , 


(2. 
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! 


where  the  volume  V is  bounded  by  the  surface  S and  n is  the  outward 
pointing  unit  surface  normal.  Q will  be  identified  with  a field 
quantity  while  P will  be  related  to  a Green's  function. 

Since  the  fields  are  vectors,  the  Green's  function  must  be  a 
dyadic  vector  operator.  The  basic  operations  and  properties  associated 
with  dyads  are  given  by  Tai  [1971,  pp.  46-48].  The  dyadic  version  of 
Green's  second  identity  is 


VxVxQ  - (VxVxP) 


0]dv 


[P  x V x Q + (V  x P)  x Q]dS  . 


(2.6) 


This  follows  from  (2.5)  by  setting  Q = Q • a where  a is  an  arbitrary 
constant  vector. 


E - Field  Formulation 

To  invert  the  wave  equations  (2.3)  the  following  dyadic  ana- 


logues of  (2.1)  and  (2.3)  are  needed: 

V x 8 (r |r')  = 5h(r|r')  (2.7a) 

V x 5h(f|f')  = ! 6(r-r')  + k2  8e(r|r')  (2.7b) 

V x V x 8 (r | r ' ) - k2  8o(r|r')  = ! 6(f-r')  (2.7c) 

0 6 

V x V x 8h(r|r')  - k2  8^(r|r')  = V x (I  5(r-r'))  (2.7d) 


a 


i jT 


\ 


Inserting  E and  Gg  into  (2.6)  yields 


E(?')  = 


f 

‘ 

. 

j(r)  • 5e(r|r')dv 


(2.8) 


| [n  x E (r) ) • Sh(r|r')  + n *(H(r)  x Se(r|r’))]dS 
S 


In  Appendix  A it  is  shown  that 


5e(r' |r)  = Se(r|r')  , 


(2.9) 


where  signifies  the  transpose  of  the  dyad.  Interchange  of  the 
primed  and  unprimed  coordinates  in  (2.8)  and  use  of  (2.9)  give 


E(i)  = 


G (r | r * ) • J(r')dV' 


(2.10) 


[(n'  x E(r'))  • Gh(r'|r)  + CG  (r|r')  x n')  • H(r')]dS' 


Since  the  tangential  electric  field  is  zero  on  the  perfectly 


conducting  surface  Sc  and  the  radiation  condition  at  infinity  holds  for 


the  fields  and  Green's  functions,  equation  (2.10)  becomes 


E(?)  = 


5 (r |r')  • j(f')dV'  - 


(S'  x E(r'))  • 5,  (r'  |?)dS', 


(2.11a) 


where 


n x 5e(r|r')  = 0 


on  S , 
c 


(2.11b) 
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n x E(r)  =0  on  Sc-A  (2.11c) 

A is  the  aperture  portion  of  Sc  and  Sc~A  is  all  of  Sc  except  the  aper- 
ture surface.  Note  that  the  fact  that  (2.9)  and  (2.11b)  imply 

S(r|r,)xn'=0  onS 

e 1 c 

has  been  used  to  obtain  (2.11a). 


H - Field  Formulation 

The  inversion  of  the  magnetic  source  wave  equations  (2.4)  is 
accomplished  similarly.  The  dyadic  analogues  of  (2.2a-b)  and  (2.4)  are 


V x gh(r|r')  = ge(r|r')  (2.12a) 

V x ge(r|r')  = k2  gh(r|r')  + f 6(r-r')  (2.12b) 

V x V x gh (r | r ' ) - k2  gh(r|r')  = f 6(r-r)  (2.12c) 

V x V x g (r|r')  - k2  g (r|r')  = V x (I  6(r-r'))  (2.12d) 

6 6 


Using  duality  between  the  electric  and  magnetic  source  driven  equations 
one  obtains 


H(?) 


> 

gh(f|r')  • M(r')dV' 
V 


f[(n’  x H(r ' ))  • L(r'|r) 
S 


(gh(r|r')  x n') 


E(r')]dS’  . 


+ 


(2.13) 
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To  obtain  the  above  one  needs 


Ih(r' I?)  = gh(r|r') 


(2.14) 


which  is  shown  in  Appendix  A.  The  radiation  condition  at  infinity  plus 


n x ge(r|r')  = 0 


on  Sc  , 


(2.15a) 


and 


n x E(r)  = 0 


on  S -A 
c 


(2.15b) 


imply 


H(?) 


gh(f|r')  • M (r ' ) dv ' - Jj  gh(r|r')  • (ii'  x E(r'))dS' 

A 

(2.15c) 


where  A is  the  aperture  portion  of  Sc< 


Dyadic  Green's  Functions  for  a Rectangular 
Cavity,  Electric  Sources 

Green's  functions  for  electric  sources  in  a rectangular  cavity 
are  constructed,  first  by  use  of  a dyadic  vector  potential  and  then  from 
the  L,  N,  M vector  wave  functions.  The  cavity  dimensions  are  a,  b,  and 
c in  the  x,  y,  and  z directions  respectively.  By  (2.7)  5^  is  divergence 
free.  Thus  it  can  be  expressed  as  the  curl  of  a vector  potential  5^. 
Maxwell's  equations  yield  [Jackson  1962] 

- (V2  + k2)  SA(r|r')  = ! 6 (r-r ' ) (2.16a) 


- (V2  + k2)  $ = - p-  V S(r-r') 


(2.16b) 


V • 5 (r|r')  = - k2  $ 


(2.16c 


where 


S^Crlr')  = V x SA(r|r')  (2.16d^ 

and 

8e(r|r')  = SA  - V $ = (!  + ^-)5a  (2.16e; 

k 

Since  the  dyadic  Green's  functions  are  field  responses  to 
sources  at  r',  they  satisfy  electromagnetic  field  boundary  conditions 
on  the  perfectly  conducting  surface  Sc,  provided  r'  is  not  on  Sc>  In 
all  considerations  of  boundary  conditions  in  this  section  it  is  assumed 
that  ?'  is  not  on  Sc.  Jones  [1974]  has  shown  that 

V • (n  x E)  = iuiun*H  (2.17) 

on  a surface  away  from  sources;  Vq  is  the  surface  divergence.  Thus 

n x S (rlr')  =0  on  S (2.18) 

e 1 c 

implies 

n • (V  x GA(r|r'))  =0  on  Sc  . (2.19) 

Note  that  (2.16e)  and  (2.18)  imply 

n x 5.(r|r')  = n x V $(r|r')  on  S . (2.20) 
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If  Dirichlet  boundary  conditions  are  imposed  on  <5,  one  obtains 


- (V2  + k2)  §A (r | r ' ) = I S(r-r') 


n x §ACr|r ' ) = 0 

> 

V • 8A(r|r')  - 0 


on  S 

c 


(2.21a) 


(2.21b) 


From  a vector  Green's  theorem  [Morse  and  Feshbach  1953,  p.  1768]  one 
may  show  that  the  operator  determined  by  (2.21)  is  self-adjoint.  Thus 
each  component  §A  • x,  GA  » ?,  and  §A  • z may  be  given  in  terms  of 
eigenfunction  expansions.  Seidel  [1977]  has  shown  that  for  r not  equal 
to  r'  the  operations  indicated  by  (2.16d)  and  2.16e)  on  §A  may  be 
performed  term  by  term  to  yield  and  Ge-  The  dyads  5 A and  Gg  are 
given  in  Table  1. 

Not  only  is  Gg  correct  away  from  the  source;  it  is 
distributionally  correct  [Stakgold  1968,  p.  11].  This  fact  can  be 
readily  demonstrated  by  reconstructing  Ge  in  terms  of  L,  M,  and  N 
vector  wave  functions  [Morse  and  Feshbach  1953,  p.  1766].  L,  M,  and 
N are  eigenfunctions  of  the  vector  Helmholtz  equation  C2.21).  Although 
(2.21)  has  a complete  set  of  eigenfunctions,  we  will  check  a complete- 
ness relation  to  insure  that  all  eigenfunctions  have  been  obtained. 

Since  M and  N are  divergence  free  they  are  also  eigenfunctions 
of  the  operator 


= V x V x E 


(2.22a) 


(cc)  (ss)(ss) 


The  cavity  has  lengths  a,  b,  and  c in  the  x,  and  £ directions  respectively. 
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with  boundary  condition 


n x E = 0 


(2.22b) 


L being  curl  free  is  in  the  null  space  of  £.  Therefore  the  vector  wave 
functions  are  also  the  eigenfunctions  of  £ and  may  be  used  to  construct 


The  vector  wave  functions  are  given  in  Table  2.  One  may  check 
that  over  the  cavity  volume  these  functions  are  mutually  orthogonal  and 
satisfy 


111  |LdS> 


2.v  _ vz  abc 

dV  ' KJlmn  e e e. 

n m Z 


(2.23a) 


HI  I*d5>l2‘v-lll 


■rr  (k  2 + k 2)abc 

tv-  m^ii’dv  — 

Z n m 


(2.23b) 


for  all  non- zero  wave  functions.  In  the  above 


1 , i = 0 


2 , i + 0 


Following  Tai  and  Rozenfeld  [1976],  we  write 


! 6 (r-r ' ) = l 

Z,  m,n 


V.£n 


* Sds>  Sds'»'  • (2-24) 

+ k 


K + *9 


i 


V (s  s sj 
x y z 


■ VxVU*  * lAVi?  * 


k s s c z 
z x y z 


= V x(cxcysz2)  = - kycxsyszx  + kxsxcys^ 
= 1_  V x V x (sxsycz2) 


k k k k k"  + 1 

-2-i.  c s s x - — — s c s y + (— 

KHmn  x y K&mn  y K«-mn 


-^-1  s s c 
1 x y z 


cxsys,  = cos(kxx)  sin(kyy)  sinCkzz) 


etc. 


k 


x 


miT 

a 


■^tT  ir  2 


T • K 


£mn 


m = 0,1,2,... 


n = 0,1,2, .. . 


and  l = 0,1,2,... 
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Inserting  the  vector  wave  functions  from  Table  2 into  (2.24)  gives 


f 6 (r-r ' ) 


£,m,n=0 


£*£m£n 
abc 


[(cc)x(ss)y(ss)z  xx 


+ (ss)x(cc)y(ss)z  yy  + (ss)x(ss)y(cc)z  22]  , (2.25) 

where  (cc)x(ss)y(ss) z 

= cos(kxx)  cos(kxx')  sin(kyy)  sinfk^y')  sin(kzz)  sin(kzz')  , 

etc. 

The  above  equation  is  known  to  be  true  by  scalar  spectral 
theory.  Note  that  dividing  each  term  of  (2.25)  by  K2^  - k2  yields 
5.  in  Table  1.  Since  L,  M,  and  N are  also  eigenfunctions  of  the 
operator  £ in  (2.22) 


8 (r|r’)  = 


” AiS* 

^ abc 


JL,m,n 


- k2  K2 


£mn 


') 


(K2  + Kz)  (K 

V x yJ  '■ 


2 

Jlmn 


k2) 


(2.26) 


Inserting  the  vector  wave  functions  yields  the  same  as  in  Table  1. 
Thus  the  vector  potential  approach  gives  the  correct  answer  even  at  the 
source  point.  The  sum  over  the  L functions  yields  a term  proportional 
to  the  longitudinal  6-function  [Morse  and  Feshbach  1953,  p.  1781]. 

The  second  term  determines  a dyad  whose  entries  are  square  integrable 


17 


over  the  cavity  volume.  This  characterization  of  the  singularities  of 
Gg  may  prove  useful  in  later  numerical  considerations. 

Dyadic  Green's  Functions  for  a Rectangular 
Cavity,  Magnetic  Sources 

Dyadic  Green's  functions  for  magnetic  sources  in  a rectangular 
cavity  are  now  constructed.  The  cavity  dimensions  are  the  same  as  in 
the  last  section.  By  (2.12)  gg  is  divergence  free  and  may  be  expressed 
as  the  curl  of  an  electric  vector  potential.  As  before  Maxwell's 
equations  yield 


- (V2  + k2)  g_(r|r')  = ! «(r-r') 


(2.27a) 


(V2  + k2)  $(r|r')  = V6lr'r.^ 


(2.27b) 


V • g_(r|r')  = - k2  $ 


(2.27c) 


where 


ge(r|r')  = V x g_(r|r') 


(2 . 27d) 


gh(r|r')  = gp(r|f ' ) - V $(r|r')  = (!  + ^j-)  gp(r|r')  . (2.27e) 


For  the  6- source  away  from  the  perfectly  conducting  surface  we  have 


n x (V  x g_)  = 0 


on  S 

c 


(2.28) 


Thus,  (2. 27e) , and  (2.17)  yield 


n * h - 37T  = n * *h  = 0 on  Sc 


(2.29) 


If  homogeneous  Neumann  boundary  conditions  are  imposed  on  $ one  obtains 


- (V2  + k2)  IF(r|r')  = ! 6(r-r') 


(2.30a) 


n x (V  x gF(r|r'))  = 0 


n • g (r|r’)  = 0 


on  S 

c 


(2.30b) 


Using  a vector  Green's  theorem  [Morse  and  Feshbach  1953,  p.  1768]  one 
may  show  that  a self-adjoint  operator  is  associated  with  (2.30).  Thus 
gc  may  be  expanded  in  terms  of  a complete  set  of  eigenfunctions.  This 
result  is  given  in  Table  3.  Again  Seidel  [1977,  p.  12]  has  shown  that 
gg  and  g^  may  be  obtained  by  differentiating  gp  term  by  term  provided  r 
is  not  equal  to  r*. 

We  now  show  that  a complete  set  of  eigenfunctions  of  (2.30)  are 
the  vector  wave  functions  L,  M,  and  N.  The  vector  wave  functions  are 
given  in  Table  4.  Since  M and  N are  divergence  free,  they  are  also 
eigenfunctions  of  the  operator 


9H  = V x V x H 


(2.51a) 


with  boundary  condition 


n x (V  x H)  = 0 


(2.31b) 


Table  3.  Magnetic  Source  Dyadic  Green's  Functions 


Table  4.  Vector  Wave  Functions  for  Magnetic  Sources* 


Eimn  = 7(cxcycz) 


= " kxsxcyCz*  * kyVyV  " kzCxCyV 


fllmn  = 7 x (sxSycz2) 


= ksccx  - k c s c y 
yxyz  xxyz7 


= v V x V x (c  c s z) 

*’mn  Kimn  x y 2 


k k k k k2  + k2 

v svcvC7^  " v Z c s c y + (_— -)  c c s z 

Knm  x y z ^ * y hm  x y 2 
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The  curl  free  L vectors  are  in  the  null  space  of  9.  It  may  be 
checked  that  the  L,  M,  and  N vectors  are  mutually  orthogonal  and 
satisfy 


k2  abc 
2-mn  e £ £„ 

m n Jo 


(2.32a) 


(k2  + 
x 


k2) 


abc 
£ £ E. 

m n SL 


(2.32b) 


for  all  non-zero  vectors.  Following  Tai  and  Rozenfeld  [1976]  the 
completeness  relation  may  be  expressed  as 

KJmn 

, W ^ + W 

k2  + k2 
x y 

II 

This  may  be  verified  by  inserting  the  vector  wave  functions  from 
Table  4. 

(ss)x(cc)y(cc)z  0 0 

0 (cc)x(ss)y(cc)z  0 

„ (cc)  (cc)  (ss) 

0 0 ^ Jx  JyK  ‘ z 


w 

I 5(r-r')  = I e0E  e 

abc„  u n SL  m n 
£.,m,n=0 


(2.33) 


Ifi(r-r') 


- A Z £ £ £ 

abc  Jl,m,n  X,  n m 


(2.34) 


Since  we  have  shown  the  vector  wave  functions  to  be  eigenfunctions  of 
the  operator  9 in  (2.31), 
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I (r|r')  .-i-Jc.ee  1 

1 abc  l m n ,_2  „2 


k2  K? 


(2.35) 


The  first  term  is  proportional  to  the  longitudinal  dyadic  delta  func- 
tion, while  the  second  part  determines  square  integrable  functions  over 
the  cavity  volume.  Substitution  of  the  vector  wave  functions  from 
Table  4 in  (2.35)  yields  g^  in  Table  3.  Thus  the  expression  for  g^  is 
also  valid  at  r equals  r'. 


The  Longitudinal  Component  of  the 
Electric  Green's  Dyad 

This  section  is  concerned  with  two  intimately  related  topics, 
the  characterization  of  the  electric  Green's  dyad's  singularity  and  the 
effect  of  neglecting  the  longitudinal  vector  wave  functions  in  the 
dyad's  eigenfunction  expansion. 

Hansen  [1935]  first  used  spherical  coordinate  vector  wave 
functions  to  expand  a free  space  potential  Green's  function.  The 
expansion  given  was  an  eigenfunction  expansion  in  the  angular  variables 
and  closed  form  in  the  radial  variable.  Hansen  found  that  the  corre- 
sponding electric  field  expansion,  exterior  to  a finite  sphere  contain- 
ing all  sources,  contained  no  longitudinal  vector  wave  functions. 
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Morse  and  Feshbach  [1953]  use  vector  wave  functions  to  represent 
a potential  dyadic  Green's  function  as  a full  eigenfunction  expansion. 
The  unit  delta  dyadic  is  broken  into  transverse  and  longitudinal 
components.  The  transverse  component  is  divergence  free  and  the  longi- 
tudinal component  is  curl  free.  It  is  implied  [Morse  and  Feshbach  1953, 
p.  1781]  that  away  from  the  source  the  transverse  and  longitudinal 
components  of  the  unit  delta  dyadic  vanish  identically. 

Tai  [1971]  has  no  longitudinal  vector  wave  functions  in  the  full 
eigenfunction  expansion  of  the  electric  Green's  dyad.  Subsequently  he 
concludes  [Tai  1973]  that  there  is  a missing  longitudinal  vector  wave 
contribution,  specifically,  - zz  5(r-r')/k2. 

Collin  [1973]  takes  a different  approach,  correcting  for  the 
incompleteness  of  the  E and  H model  expansions  in  a waveguide.  His 
correction  term  is  the  same  as  Tai's  zz  delta  function  term.  In  his 
conclusion  Collin  alludes  to  a private  communication  with  Tai  and 
implies  that  this  2z  delta  correction  term  must  be  added  to  the  trans- 
verse vector  wave  functions  to  obtain  the  complete  dyadic  Green's 
function. 

Howard  [1974]  uses  a cylindrical  vector  wave  function  expansion 
(the  sin  m4>  terms  should  be  added  for  completeness)  to  produce  a closed 
form  expression  for  the  longitudinal  vector  wave  function  contribution 
to  the  full  Green's  dyadic  eigenfunction  expansion.  This  term  does  not 
identically  vanish  away  from  the  source,  apparently  contradicting  the 
previously  cited  literature. 

Tai  and  Rozenfeld  [1976]  correctly  state  that  the  longitudinal 
vector  wave  functions  in  the  full  dyadic  expansion  contribute  not  only 
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at  the  source  point  but  away  from  the  source  as  well.  But  then  later 
they  seem  to  imply,  as  Tai  [1973],  that  the  longitudinal  vector  wave 
function  contribution  is  a - zz  6(r-r')/k2  term.  These  two  statements 
are  in  contradiction. 

Several  papers  deal  with  the  singular  nature  of  the  electric 
Green's  dyadic.  Van  Bladel  [1961]  expresses  the  free  space  electric 
Green's  dyadic  in  terms  of  a spherical  principal  value  plus  a delta 
function  type  singularity.  Rahmat-Samii  [1975]  extracts  different 
singularities  from  cavity  and  waveguide  dyadics.  Yaghjian  [1978]  used 
"principal  volumes"  and  delta  function  type  singularities  to  character- 
ize the  electric  Green's  dyadic's  singularity.  He  notes  that  different 
"principal  volumes"  and  their  associated  different  delta  type 

singularities  have  a sum  total  singularity  which  is  the  same  irrespec- 
tive of  the  "principal  volume"  chosen. 

Some  of  the  confusion  in  the  literature  may  be  rectified  by 
clearly  stating  which  longitudinal  vector  wave  functions  are  being  used. 
However,  the  papers  of  Howard  [1974],  Tai  [1973],  and  Tai  and 
Rozenfeld  [1973]  still  remain  in  contradiction.  Second  it  is  physically 
surprising  that  Rahmat-Samii  [1975]  has  extracted  different  singulari- 
ties for  cavity  and  waveguide  problems. 

It  is  rigorously  shown  that  the  dominant  singularity  of  the 
electric  Green's  dyad  is  proportional  to  the  longitudinal  component  of 
the  unit  delta  dyad.  Contrary  to  Morse  and  Feshbach  [1953],  this  term 
does  not  vanish  identically  away  from  the  source.  The  longitudinal 
vector  wave  functions  in  the  full  eigenfunction  expansion  of  the  Green's 
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dyad  not  only  contribute  at  source  points  but  away  from  the  source  as 
well.  The  full  eigenfunction  expansion  is  useful  for  two  reasons. 

First  spectral  theory  will  guarantee  the  completeness  of  the  eigen- 
function expansion.  Second  if  one  desires  to  compute  the  Green's  dyad 
at  different  spacial  locations,  the  full  expansion  allows  one  the 
flexibility  to  close  in  a manner  insuring  the  greatest  exponential  con- 
vergence. Howard's  central  result  for  free  space  [Howard  1974,  eq.  13] 
will  be  mathematically  verified  and  shown  to  be  unique.  The  apparent 
controversy  on  the  effect  of  the  longitudinal  vector  wave  function 
contribution  will  be  resolved. 

The  free  space  electric  Green's  dyadic  is  studied  first. 

Results  are  later  generalized  to  cavity  and  waveguide  problems.  The 
free  space  electric  Green's  dyad,  Gg,  satisfies  [2. 7c)  plus  the 

Sommerfeld  radiation  condition.  G is  obtained  by  vector  wave  function 

e 

expansions.  Expansion  of  the  unit  delta  dyadic  yields  [Morse  and 
Feshbach  1953] 


over  a continuous  or  discrete  set  of  m.  The  first  term  is  the 
longitudinal  component  of  the  delta  dyad;  the  second  is  transverse.  By 
spectral  theory 


r 


- D (r-r') 

G frir')  = 

e k2 


M (?)  M*  (r ' ) N (r)  N*(r')  . 

r m m^  ' mk  y nr  1 

^ A'  + A"  ,2  v2 

m m k - k 

m 


C2.37) 


One  may  show  that  the  second  term  is  square  integrable  over 
finite  volumes  by  noting  the  corresponding  cavity  sum  is  square  inte- 
grable and  that  the  difference  of  the  cavity  and  free  space  sum  is 
analytic.  Thus  the  dominant  singularity  of  Gg  is  proportional  to 
DL(r-r'),  the  longitudinal  component  of  the  unit  delta  dyadic.  Morse 
and  Feshbach  [1953]  indicate  that  D^(r-r')  is  identically  zero  for  r not 
equal  to  r'.  This  implies  that  the  L functions  may  be  neglected  away 
from  the  source.  However  we  proceed  to  show  that  D^(r-r')  does  not 
vanish  identically  away  from  the  source.  This  is  done  by  two  indepen- 
dent proofs. 

First  a proof  by  Helmholtz's  theorem  is  given.  Following 
Howard  [1974] 


I fi(r-r')  = V2  i——  = V x V x TV  — (2.38) 

4tt  I r-r ' I 4ir|r-r'|  4tt  I r-r  * I 


The  first  term  is  transverse;  the  second  longitudinal.  One 


expects  that 


D^(r-r')  = - W 


4iT|r-r' 


(2.39) 
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A mathematical  proof  of  (2.39)  is  now  given.  It  is  also  shown  that  this 
dyadic  distribution  is  unique. 

First  note  that  any  dyadic  distribution  is  uniquely  determined 
by  its  operation  on  infinitely  differentiable  testing  vectors  F(r)  with 
compact  support.  This  may  be  seen  by  putting  the  dyad  in  matrix  form 
and  F in  vector  form  and  the  use  of  scalar  distribution  theory.  The 
Helmholtz  theorem  (Morse  and  Feshbach  1953]  shows 

F (r)  = V $ + 7 x A ; V • A = 0 (2.40a) 

where 


*(r)  « - 


(f  V • F(r') 


4tt  r-r' 


dv' 


and 


(2.40b) 


rrr 

Mr)  = 


V-2LjSZll.  dv. 

4ir|r-r'  | 


(2.40c) 


Note  that  surface  terms  in  (2.40a)  and  (2.40b)  are  not  present  since  F 
has  compact  support.  Thus  the  longitudinal  part  of  the  test  vector  F 
is 


fl(?)  = - V 


ZLlMH  dv. 

47r|r-r'  | 


(2.41a) 


Kellogg  [1953]  shows  the  gradient  may  be  brought  inside  the  integral; 


fi 


thus 


T 


^ L 


fl(?)  = - 


JJ 


[V*  • F(r')]  V 


dv' 
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(2.41b) 


4tt  r-r' 


From  distribution  theory 


FLW 


fl 


F(r')  • W 


dv' 


(2.41c) 


4ir  r-r' 


Since  a distribution  is  determined  uniquely  by  its  operation  on  testing 
vectors  F(r), 


DL(r-r')  = - W 


1 


(2.42) 


4ir  r-r ' 


is  uniquely  determined.  Direct  computation  shows  that  D^(r-r')  does  not 
vanish  identically  away  from  the  source. 

An  alternate  constructive  proof  of  the  above  result  (2.42)  is 
now  given.  Expansion  of  the  free  space  delta  dyadic  in  terms  of 
cartesian  coordinate  vector  wave  functions,  (2.36),  and  orthogonality 
yields 


f 5 (r-r ' ) = B^(r-r')  + 5^,(r-r') 


(2.43a) 


where 


B (r-r')  = 

*-*  /"w*  3 


rr 


(2tt) 
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dk  dk  dk 
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(2.43b) 
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x y z 

(2.43c) 
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(2.43d) 


and 


4;  = e~*K  ^ , K = k x + k y + k z 
x y z 


By  direct  computation 
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k2  + k2  + k2  x ^ 

x y z 


(2.44a) 


(2.44b) 


= - W 


4tt  I r-r ' 


(2.44c) 


which  is  identical  to  (2.39). 

Although  the  previous  discussion  was  restricted  to  free  space 
one  may  suspect  analogous  results  to  hold  for  cavity  and  waveguide 
problems.  This  is  indeed  the  case  as  is  now  shown.  The  electric  field 
cavity  dyadic  is  given  by  (2.16e)  where  5^  satisfies  (2.21).  The 
superscript  "c"  has  been  added  to  cavity  dyads  to  avoid  confusion  with 
free  space  dyads.  The  boundary  value  problem  (2.21)  is  self-adjoint 
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and  has  a complete  set  of  associated  eigenfunctions.  These  same  wave 
functions  may  be  used  to  expand  the  delta  dyad  and  as  in  (2.36)  and 
(2.37).  Although  Rahmat-Samii  [1975]  extracted  a delta  singularity  from 
Sg,  its  true  dominant  singularity  is  B^(r-r')  defined  analogously  to 


(2.36). 


As  a device  to  study  D^(r-r')  we  consider  5^  in  (2.16e)  with  k 


set  equal  to  zero,  5^  . Expanding  in  vector  wave  functions 


(r  r’) 


« " A;  k»  * A-  k» 


(2.45) 


From  (2.16e) 


v2  fe(?i?')  - A—  | = ° 

\ 4u|r-r' | J 


(2.46) 


interior  to  the  cavity  or  waveguide.  From  the  ellipticity  of  the  above 
operator  it  follows  that  the  difference  of  the  two  potentials  is 
analytic  [Treves,  1975].  However 


W • §£(^')  - 


4ir  r-r' 


Lm^  Lm^’ 


Dj (?-?') 


(2.47a) 


= D^(r-r ' ) - Dl  (?-?') 


(2.47b) 
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To  obtain  (2.47a)  term  by  term  differentiation  was  used  which  is 
distributional ly  correct  [Stakgold  1968].  Since  the  left  side  of 
(2.47a)  is  analytic  it  follows  that  D^(r-r')  has  the  same  singularity 
as  D^fr-r1).  Thus  D^(r-r’)  does  not  vanish  identically  for  r/r'. 

Although  the  above  results  are  for  r,r'  in  the  cavity  or  waveguide 
interior,  by  use  of  an  image  term  one  can  obtain  similar  results  for 
r'  on  a wall. 

Different  representations  of  the  electric  Green's  dyadic  are 
now  considered.  The  electric  Green's  dyad  may  be  expressed  as  an  eigen- 
function expansion  in  two  variables  and  closed  form  in  the  third.  As 
previously  stated,  this  is  the  form  in  which  Hansen  [1935]  first  intro- 
duced the  vector  wave  function  expansions  and  it  was  found  that  no 
longitudinal  vector  wave  functions  were  needed  outside  of  a sphere 
containing  all  sources.  It  will  now  be  shown  that  the  effect  of 
neglecting  the  longitudinal  wave  functions  in  this  form  is  precisely 
the  zz  delta  function  correction  term  [Tai,  1973;  Tai  and  Rozenfeld  1976; 
Collin,  1973].  Although  free  space  is  considered,  results  for  wave- 
guides follow  analogously. 

Closing  (2.44b)  in  the  complex  kz  plane  and  then  differentiating 
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Closing  the  transverse  part  of  5g  (2.37)  in  the  complex  plane 
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where  k^  = A2  - k2  . Recalling  (2.37),  we  see  that  the  second  term  of 
(2.49),  the  static  pole  contribution,  cancels  with  the  contribution  from 
the  first  term  of  (2.48).  The  remaining  terms  of  Gg  are  just  the  closed 
form  in  z,  eigenfunction  in  x and  y,  M and  M expantion  function  contri- 
butions plus  the  famous  - 22  6(r-r')/k2  term. 

In  conclusion  the  true  dominant  singularity  of  the  electric 
Green's  dyadic  is  proportional  to  the  longitudinal  component  of  the  unit 
delta  dyadic.  Contrary  to  Morse  and  Feshbach  [1953],  this  term  does  not 
identically  vanish  away  from  the  source.  Rahmat-Samii  [1975]  only 
extracted  part  of  the  singularity  from  his  cavity  and  waveguide  dyads , 
their  true  dominant  singularities  are  the  same  away  from  the  cavity  or 
waveguide  walls.  At  the  walls,  image  terms  must  be  considered.  The 
longitudinal  L vector  wave  functions  in  the  full  eigenfunction  expansion 
of  the  electric  Green's  dyad  contribute  not  only  at  the  source  but  away 
from  the  source  point  as  well.  The  full  eigenfunction  expansion  is  use- 
ful for  two  reasons,  one  theoretical  and  one  practical.  First  spectral 
theory  guarantees  the  completeness  of  the  electric  Green's  dyadic 
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expansion.  Second  if  one  desires  to  compute  the  Green's  dyad  away  from 
the  source,  one  has  the  flexibility  to  close  in  the  variable  which  gives 
fastest  exponential  convergence. 

The  confusion  regarding  the  longitudinal  vector  wave  function 
contribution  in  the  literature  has  been  resolved.  In  Tai's  early  formu- 
lation [Tai,  1971]  the  missing  L function  contribution  is  non- zero  away 
from  the  source.  Only  when  the  Green's  function  is  displayed  as  a two- 
dimensional  eigenfunction  expansion  with  closed  form  in  the  third 
dimension  is  the  L function  contribution  completely  given  by 
- zz  S(r-r')/k2.  It  is  crucial  to  note  that  if  this  partially  closed 
representation  is  used,  the  22  delta  function  term  is  not  the  total 
singularity  since  the  ft  and  N function  contribution  is  highly  singular. 
The  total  dominant  singularity  is  still  proportional  to  the  longitudinal 
component  of  the  unit  delta  dyadic.  Finally  similar  techniques  may  be 
applied  to  show  that  the  cavity  dyad  (2.35)  has  the  same  type  of 
singularity. 


CHAPTER  3 


A PRELIMINARY  PROBLEM 


In  an  effort  to  gain  insight  into  some  of  the  difficulties  that 
may  occur  in  the  three-dimensional  vector  problem  of  Chapter  1,  a scalar 
two-dimensional  diffraction  problem  is  considered.  This  problem  illus- 
trates how  to  deal  with  the  field  and  kernel  signularities  involved  in 
solution  of  the  aperture  integral  equations.  As  expected,  the  dyadic 
formulation  for  this  problem  reduces  to  the  scalar  formulation. 

The  problem  considered  is  the  diffraction  of  steady  state 
electromagnetic  waves  through  an  aperture  in  a thin,  perfectly  conduct- 
ing, infinite  sheet  (Figure  2).  As  illustrated,  the  geometry  and  the 
sources  are  y independent;  thus  the  fields  also  have  no  y variation. 
Maxwell's  equations  for  electric  sources  (2.1)  decouple,  viz: 


- (V2  ♦ k2)  Ey  = Jy 


/N 


(3.1a) 


(3.1b) 


A 


BE 


(3.1c) 


The  magnetic  source  analogues  (2.2)  becomes 


- (Vz  ♦ kz)  Hy  = My 
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(3.2a) 


u 


3H 

E I 

x 3z 


(3.2b) 


3H 

E = 
z 3x 


(3.2c) 


Equations  (3.1)  and  (3.2)  have  been  called  the  "hard"  and  "easy" 
polarizations  respectively. 


The  Easy  Polarization 

An  integral  equation  for  the  tangential  electric  field  across 
the  aperture  is  developed.  Once  this  tangential  electric  field  is 
known  Green's  second  identity  [Kellogg  1953,  p.  215]  may  be  used  to 
obtain  at  any  point.  The  free  space  Green's  function  for  (3.2a)  is 
given,  apart  from  a complex  conjugate  needed  to  account  for  an  ela)T 
time  dependence,  by  Morse  and  Feshbach  [1953,  p.  811]  as 


go(x,x’,z,z')  = jr-  H^2-1  (k/(x-x')2  + (z-z')2  ) . 


Thus  the  half-space  Green's  function  satisfying  the  appropriate 
magnetic  boundary  condition  is  given  by 


g(x,x',z,z')  = ^ [H<’>  (k/(x-x’)2  + (z-z’)2  ) 


+ Hj2J  (k/(x-x')z  + (z+z')2  ) ] • 


(3.3) 


(3.4) 


Using  Green's  second  identity,  the  radiation  condition  at  infinity,  and 
continuity  of  the  tangential  field  components  across  the  aperture,  one 


obtains 


2iH  (x) 

y v J 


L, 


hJ2)  (k|x-x'|)  Ex(x')  dx' 


(3.5a) 


-L 


The  incident  field  is  given  by 


2 Hlnc (x) 

y 


g(x,x',o,z')  My(x',z')  dx'dz ' , 


(3.5b) 


where  is  all  of  Region  I.  The  above  equation  is  solved  numerically 

A 

by  approximating  E^(x)  by  pulses,  point  matching,  and  inverting  a matrix 
equation.  As  illustration,  numerical  results  are  given  for  a normally 
incident  plane  wave  illuminating  an  aperture  of  total  width  . 025A 
(Figures  3 and  4).  The  dashed  lines  give  computed  values  of  when 
expanded  in  20  pulses;  the  dotted  lines  give  E^  using  a 40  pulse 
expansion. 


The  Hard  Polarization 

Techniques  similar  to-  those  used  to  solve  the  easy  polarisation 
may  be  applied  to  (3.1)  to  obtain 


2iHinc(x) 
x v J 


= lim 
z-+o 


-L 


3z 


-HC2)  (kt/(x-x' ) 2 

2 O 


+ z2  ) 


Ey(x')  dX( 


(3.6) 


where  Hx  (x)  is  the  x component  of  the  incident  magnetic  field  in  the 
aperture.  As  z approaches  zero  the  kernel  becomes  highly  singular  and 
the  limit  and  integration  cannot  be  interchanged.  Noting  that 


i nc  • 

Figure  4.  Real  E , Easy  Polarization 
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(il_  + ll_  + k2)  H(2^  (k/(x-x')2  + z2  ) = 0 , z t 0 (3.7) 

Bx2  3z2  0 


(3.6)  becomes 


2iH^nC(x) 


lim(ii-  ♦ k2) 

**°  3x2 


hJ2)  (k/(x-x* ) 2 + z2  )Ey(x’)dx' 


(3.8) 


Johnson  and  Dudley  [1978 ] have  studied  the  aperture  electric 
field's  singularities  and  have  obtained  a method  of  numerical  solution 
of  (3.8)  with  an  interpolary  feature  that  may  be  useful  in  the  higher 
dimensional  vector  problem.  Alternately  (3.8)  may  be  solved  by  pulse 
expansion  and  piecewise  sine  testing  [Wilton  and  Butler  1976],  A 
comparison  of  Johnson  and  Dudley's  Green  function  moment  method  (GFM) 
and  Butler's  piecewise  sine  method  (PS)  is  given  for  normally  incident 
plane  wave  excitation.  The  real  and  imaginary  parts  of  the  aperture 
electric  field  versus  the  aperture  coordinate  are  given  for  total 
aperture  widths  of  .01A  (Figure  S)  and  1.8A  (Figure  6). 

Dyadic  Formulation  for  the  Easy  Polarization 
The  dyadic  formulation  for  the  easy  polarization  is  now  given 
and  shown  to  reduce  to  the  scalar  problem.  By  image  theory  the  electric 
vector  potential  dyad  (2.30)  for  the  half-space  one  or  two  (Figure  2)  is 
given  by 

gp(rjr')  * 01/  + + >/)yy  * C'P  ~ ^i)2z  , (3.9a) 


where 


-ikR1 


r = [(x-x')2  ♦ (y-y')2  + (z-z')2]1/2 


r1  = [(x-x')2  + (y-y')2  + (z+z')2]1/2 


From  (2. 27e)  one  obtains, 


, (r | r ' ) = (!  + -j1)  IF(r|r') 


Equivalently, 


g.  (r|r')  = (!  - — VV')0j>  + ^)  - 2 ip1  11 


Equation  (2.15c)  yields 


Hi(r)  = gh(f|r')  * 9 My(x’,z')dv 


gh  (r|r')  • (2  x E1(?'))ds' 


(3.9b) 


(3.9c) 


(3.9d) 


(3.9e) 


(3.10) 


(3.11) 


(3.12) 
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for  ? in  Region  1,  where  A is  the  aperture  surface.  In  Region  2, 


h2(?) 


= | gh(r|r')  • (z  x E2(r'))ds' 


(3.13) 


Noting  the  relationship  between  the  two  and  three  dimensional 
Green's  functions  [Morse  and  Feshbach  1953,  p.  1323]. 


(k*6c2+  z2) 


-ik/x2  + y2  + z2 


_oo  4ir/x2  + y2  + z2 


(3.14) 


the  volume  integral  in  (3.12)  becomes 


9 ||  My(x\z’)  jr-  [H^2)(k/(x-x’)2  ♦ (z-z')2) 


♦ H^2  (k/(x-x')2  + (z+z' ) 2) ]dx'dy' 


(3.15) 


Continuity  of  the  tangential  field  components  across  the  aperture  and  y 
independence  of  the  fields  yield 


9 ||  My(x',z’)  jt-  [H^2)(k/(x-x')2  ♦ (z-z')2) 


+ (k/(x-x')2  + (z+z')2)]dx'dz' 
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fL 


= y 


(k|x-x'  |) 

E (x’)  — ^ dx' 

x i 


-L 


x lim  (—  + k2) 
z-*o  3x2 


L H^2^  (k/(x-x')Z  + z2> 


E (X1) 

y 


dx' . 


ik2 


-L 


(3.16) 


Since  only  the  x component  of  the  electric  field  is  excited  E^  is  set 
to  zero  and  (3.16)  reduces  to  its  scalar  counterpart  (3.5). 


Dyadic  Formulation  for  the  Hard  Polarization 
The  dyadic  formulation  for  the  hard  polarization  also  reduces  to 
its  corresponding  scalar  equation  (3.8).  By  image  theory  the  magnetic 
vector  potential  dyad  (2.21)  for  the  halfspace  is  given  by 

gA(r|r')  = (ip  - \pi)xx  + 0J>  - 4'1)yy  + + ^)z2  , (3.17) 


where  ijj  and  aTe  defined  in  (3.9).  Equation  (2.16e)  may  be  manipu- 
lated into  the  form 


E^r) 


5e(*l?’) 


y j (x' ,z')dv' 


(z  x E(r')) 

A 


• V’  x Se(r' |r)ds' 


(3.18a) 


in  Region  1.  Proceeding  as  before  (3.12  - 3.15),  (3.18a)  may  be 
expressed  as 


I<CU 


r8 


. zzz  'zzrxr*"-  -yr^iasn.r'-xaz.^zss''  - ■ 


rr 


(r)  = y 


Jy(x',z') 
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[Ho(k/(x-x')  + (z-z ' ) 2)-Ho(k/(x-x' ) 2+(z+z ' ) 2)  ] dx'dz' 

_ 


(2  x ILfr'))  • V'  x 5 (r' |r)ds' 

X 6 


(3.18b) 


In  Region  2 (2.11)  yields 


E2(r)  = ||  (2  x E2(?'))  • V x 5e(i'|r)ds' 


Using  (3.14)  and  the  y independence  of  the  fields 


(3.19) 


Hlt^  = 


= x 


» H*-2-*  (k/(x-x' ) 2 + z'  2)  dx'dz' 

^y(x',Z,)  dT'  2x 


lim 
z->o  - 


[X 


az< 


E (x') 

y 


H^(k/(x-x')2  + z2)  dx' 

2l 


-L 


A 1 2 

+ y k 


L A H(2)(k/(x-x')2  + z2)  dx' 


Ex^  2i 


] (3.20a) 


Where  H (r  ) is  the  tangential  H component  as  r approaches  the  aperture 
it  3. 

point  r from  Region  1.  Similarly 

cL 


-L 


Equating  the  x components  of  (3.20a)  and  (3.20b)  yields  the  scalar 
equation  (3.6).  Since  the  x component  of  the  electric  field  is  not 
excited  it  is  set  to  zero.  One  may  check  (3.18)  or  (3.19)  to  see  that 
indeed  only  the  y component  of  i is  excited.  Thus  the  dyadic  problem 
has  reduced  to  the  scalar  problem  providing  a check  on  the  dyadic 
formulation. 

is 


^ . u 


CHAPTER  4 


INTEGRAL  EQUATIONS  FOR  APERTURE  FIELDS 
AND  OBSTACLE  CURRENTS 


The  three  dimensional  exterior-interior  coupling  problem  of 


Chapter  1 is  now  considered.  The  dyadic  formalism  of  Chapter  2 is  used 


to  develop  a set  of  coupled  integral  equations  for  unknown  aperture 


fields  and  obstacle  currents.  The  techniques  of  Wilton  and 


Dunaway  [1975]  for  aperture  penetration  through  a perfectly  conducting 


infinite  screen  are  followed. 


Application  of  (2.10)  to  Region  1 of  Figure  1 yields 


EW  = ESC(r)  + 


[z  x E(r')]  • S£S(r' |r)  ds' 


(4.1a) 


where 


esc 

E (r)  = 


5^S(r  r')  • J(r')  dv' 


Region  1 


is  interpreted  distributionally  and  the  superscript  "hs"  denotes  a 


half  space  dyad.  The  short  circuited  electric  field  Es  is  the  field 


that  would  be  present  if  the  aperture  were  not  present.  Similarly  in 


Region  2 (Figure  1)  one  obtains 


E(?)  = - j|  [z  x E(r' ) ] • 5£(?'|?)  ds'  - |J  S^(?|f) 


Obstacle 


• (n  x H(r')  ds' 


(4.1b) 
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where  n is  the  unit  normal  pointing  outward  from  Region  2 and  inward  to 
the  obstacle.  The  superscript  "c"  denotes  a cavity  dyad.  On  a perfectly 
conducting  surface  the  surface  current  equals  fi  x n [Harrington  1961, 
p.  106] . Thus 


E(r)  = - 


[z  x E(r')]  • 5h(r'|r)ds'  + 


Ge(r|r') 


Obstacle 


• J(r')  ds' 

in  Region  2. 

Application  of  the  result 

[V'  x 5 (r'|r)]~  = V x g,  (r|r')  (Appendix  A] 


(4.1c) 


(4.2) 


to  (4.1a)  plus  (2.1a)  yields 


H(r)  = HSC (r)  + V x 


JJ 


V x g^S (r | r ' ) • (z  x E(r'))  ds'  (4.3a) 


interior  to  Region  1,  where  Hsc(r)  = V x Esc(r).  Similarly  (4.1c)  and 
(4.2)  yield 


H(r)  = - V x 


V x g^(rjr')  • [z  x E(r')]ds' 


+ V x 


II 


G (r|r')  • J(r')ds' 


(4.3b) 


Obstacle 


interior  to  Region  2. 


— ■ip.«m.ni'.i.u  impij 


Now  continuity  of  the  tangential  components  of  the  aperture 
fields  will  be  invoked  to  obtain  a coupled  set  of  integral  equations 
for  the  tangential  electric  field  in  the  aperture.  An  integral  equation 
will  also  be  enforced  for  the  obstacle  surface  currents  which  are 
coupled  to  the  tangential  components  of  the  aperture  electric  field. 
Equation  set  (4.3)  and  (2.12c),  plus  continuity  of  tangential  field 
components  yield 


Ht(r)  = [lira  - k2  [(  g£S(r| 
1 r*-c+  J{  h 


I?’)  * (Ex(f')?  - Ey(r')x)ds' 


+ lim  -k' 
z->-c" 


II  ‘£(f| 


r')  * (Ex(r')y  - Ey(r')x)ds' 


5^(r|r')  • J (r ' ) ds ' ] t , 


(4.4.) 


Obstacle 


in  the  aperture  where  the  subscript  "t"  denotes  the  tangential  component. 
Application  of  (2.16d)  and  (2.27 e)  to  (4.4)  yields 


usc . . 3 3 

Hx  (r)  ■ lim  - j;  ^ 


EKCr')ds' 


+ lim  ( + k2) 

z-*c  3x2 


[-hs,xX(-|-,)  + gC.xx^j-,^  Ey(r ' )ds  ' 


Obstacle 


, ^ "25T, 


in  the  aperture  interior.  Similarly  the  result 


HSvCM  = lim  - (— + k2)  ff  [gis,yy 


z-k:  3y 


+ gp,yX (?!?')]  Hx(f')ds'  + tgpS,XX(r|r' ) 

+ gp’XX(?^')]  Ey(?')ds-  + |j  [g|G^XX(?|f')  JXC^M 


Obstacle 


a!  =A,IZ«lf) 


(4.5b) 


may  be  obtained  in  the  aperture  interior.  Equation  set  (4-5)  will  be 
decoupled  by  use  of  techniques  put  forth  by  Dunaway  [1974]. 

Since  E and  H satisfy  Maxwell's  equations  they  are  not  indepen- 
dent. In  fact  continuity  of  tangential  aperture  fields  implies  that  the 
normal  component  of  E is  continuous  across  the  aperture  [Jones  1974] . 
From  equation  set  (4.1), 

Esc(?)  = 2 • [-  lim  [[  V x g[|s (r | r ' ) • [E  (r ' )y  - E (r' )x]ds ' 

2 ^ y 


-lim  Vx  g£(f|?')  * [Ev(r ’)y  - E (r ' )£]ds ' 


Gg(r|r')  • J(r ' )ds ' ] 


(4 . 6a) 


Obstacle 


Equation  (4.6a)  may  be  simplified  through  the  use  of  (2.16e)  and  (2.27e) 


viz: 


z-*c 


[g£S,yy(r|?')  - gp,yy (?!?')]  Ex(r')ds* 


. . 3 

- lim  97 

z-*c 


[gJS'XX(?|?')  > gj;’XX(?|?’)]  Ey(i')ds' 


Jr  f 


Obstacle 


+ ^7  G^^Crli')  J Cr')  + (~~r  + *2) 


3z3y  A 


3z‘ 


• g£’ZZ(?|?')  J8(r')]  ds* 


(4.6b) 


Now  (4.5)  is  brought  into  a more  convenient  form  by  use  of 
(4.6b).  The  E term  in  (4.5a)  is  replaced  by  the  appropriate  counter- 
part from  (4.6b).  This  yields 


<v‘t 


k2) 


II 


[g5s,xxc?l?') 


♦ gcc’xx(?l?')] 


A ^ — 


E (r')ds' 
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_3 

3y 


(?) 


v2 


-C,XX,-|- 


3x3y3z  A 


(?!?’) 


J (r‘) 
x 


Obstacle 


+ — GA  (f  V*'3  + (£"T  + ^ 37 
3y23z  A y 3z2  3y 


• G^’ZZ(r|r')  J2(r')]  ds' 


- .]  I17  °A,ZZ(flf,)  J2<f,)  - al  GI,yy  V?')ld5'  ’ 

Obstacle 

(4.7a) 

In  (4.7a)  the  "V2"  operator  equals  - — + - — . Similarly  (4.5b)  becomes 

3x2  3y2 


(V2  + k2)  ||  [gjs,yy(?|?')  + gpS,yy(r|f’)]  §x(r')ds' 
A 


o Z 

• gJ’Z2(?|*-)  Jz(r')]ds'  ♦ ||  ^3l  G^’xx(r|r')  Jx(r') 

Obstacle 

" 37Ga’ZZ(?I?,)  Jz^')]ds'  • (4 • 7b) 


Equation  set  (4.7)  is  desirable  not  only  for  the  partial 
decoupling  effect,  but  also  for  the  presence  of  the  harmonic  operator 
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(V*  + k2).  Although  (4.5)  plus  the  correct  obstacle  surface  currents 
is  sufficient  to  uniquely  determine  the  electromagnetic  fields,  (4.7) 
is  not.  The  aperture  fields  must  be  subject  to  the  constraint 
equation  (4.6b). 

Again  using  Dunaway's  technique  [Dunaway  1974]  the  harmonic 
operators  of  (4.7)  will  be  inverted.  To  invert  the  harmonic  operators 
consider  the  two  dimensional  Helmholtz  problem  described  in  Figure  7. 

The  problem  is  given  as 

v(p),  pe  a 

- (V2  ♦ k2)  $(P).  = { (4.8) 

0 , pe  A2 

where  $(p)  is  assumed  continuous  across  the  boundary  and  satisfies  the 
radiation  condition  at  infinity.  Use  will  be  made  of  Green's  second 
identity  in  two  dimensions  given  by 

(ip  V2  $ - <f  V2  Mds  = <’0J>  - * !jj}d£  (4.9) 

[Harrington  1961].  Also  use  is  made  of  the  two  dimensional  Green's 
function  satisfying  the  radiation  condition  at  infinity.  This  is  given 
by 

H(2)(k|p-p'|) 

g0(p|p"')  - - n (4.10a) 

and  satisfies 
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Application  of  Green's  second  identity,  (4.9),  to  A2  of  Figure  8 yields 

| tg0(p|p')  H-,  CP')  - $(p')  g0(p|p,)]di'  = 0 (4.11a) 

3A1 

provided  p is  interior  to  A^  where  n is  the  outward  pointing  unit  normal 
from  Aj.  Similarly  application  of  (4.9)  to  Aj  yields 


•J  8o(SI' 


p')  v(p')ds'  + [g„(p|p') 


3gn(p|p') 

- *(P'’)  — ^ ]dV  . 


(4.11b) 


Under  the  assumption  that  <f(p)  is  continuous  across  the  boundary, 
subtraction  of  (4.11a)  from  (4.11b)  yields 


■ | 80(»l 


p')  v(p')ds'  + 


j 


p')  e(p')dr  (4.12) 


where 

the  subscripts  "A^"  and  "A^"  denote  limiting  values  as  p'  approaches  the 
boundary  from  the  respective  region,  and  p belongs  to  the  interior  of  A^. 
Since  the  aperture  fields  satisfy  Meixner's  edge  condition  [Van 
Bladel  1964]  , one  may  show  that  the  functions  operated  on  by  the 
Helmholtz  operators  (4.7)  satisfy  the  appropriate  continuity  conditions 


■ ilGA,yyc?|?,)  y?,^ds' 


and  ? belongs  to  the  aperture  surface  A.  Equation  (4.7b)  becomes 


[ [gpS,yytfl?’)  ♦ gp,yy C?|r')]  Extf')ds' 

A 

= g0(?|r')  \JJx(r')d£'  + g0(r|r') [-H*c(r') 


- -g|  E|cCr‘ ) + Tx(r')  + T4(f')]ds'  , 


(4.13b) 


where 


T_(r)  = — [f  [ ~~  ’ xx (r | i ' ) J (?') 

3 v2  JJ  avz»7.  A 


k2  J . 3x  3z 
Obstacle 


33  rc. 


3 , 32 


Jycf) 


• cJ[’ZZtf|?')  Jz(r')]ds'  , 


Vf>*  'll  (wGAXX^f"1  Vf'> 

Obstacle 


^IGI'“C5|5M  a*’  . 


and  r belongs  to  the  aperture  surface. 

Equations  (3.13a)  and  (5.13b)  are  coupled  by  means  of  the 
obstacle  currents  and  the  auxiliary  unknowns  ip  and  ip  . Before  proceed- 

X y 

ing  to  express  the  constraint  equation  in  terms  of  obstacle  currents  and 
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the  auxiliary  unknowns,  it  is  convenient  to  express  the  short  circuited 
aperture  fields  in  terras  of  the  incident  fields. 

Figure  8 illustrates  the  incident  plane  wave  electric  field. 

The  incident  electric  field  is  given  by 


Einc(x,y,z)  = [(Ejnc  cos0  cosi}>  - E^nc  sin<J>)x 


+ (Eq110  cos0  sin<J>  + E^n  cos<J>)y 


ginc  s^nQ  g]  sin0  cos<})  x gik  sin0  sintj)  y gik  cos0  z 


Since  the  short  circuited  fields  satisfy 


(4.14) 


V x ESC  = HSC  , 


(4.15) 


one  obtains 


usc  9 esc 

n “ c 

x 9y  z 


9 isc 
9z  y 


2 ^ EinC  | 

9z  y 1 z=c 


(4.16a) 


„sc  _9  esc  _ _9  esc  2 JL  p^nc 

y + 9x  z ~ 9z  x 9z  x | z=c 


(4.16b) 


in  the  aperture.  These  equations  may  be  used  to  replace  expressions  for 
short  circuited  aperture  fields  in  equation  set  (4.13)  by  their  incident 
field  counterparts. 
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Although  (4.5)  plus  the  boundary  condition  that  n x E = 0 on 
the  perfectly  conducting  obstacle  (4.1c)  are  sufficient  to  determine 
the  electromagnetic  fields,  the  constraint  equation  (4.6b)  must  be 
enforced  if  the  more  tractable  set  (4.13)  is  to  replace  (4.5) 

[Mittra  et  al.  1973].  The  constraint  equation  (4.6b)  will  now  be  given 
in  terms  of  the  auxiliary  unknowns  plus  the  obstacle  currents. 

Substitution  of  (4.13)  into  (4.6b)  yields 


= ‘ if  / Vf,)dA' 

3A 
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+ (—  + k2)  G?’  zz (r  | r' ) J (?')]ds< 


3zx 


Next  the  integration  by  parts  formula. 


3v  , 

u ■* — dx  = 
3x. 

1 


v dx  + | u v n 6.V 

3xi  l xi 


where  n is  the  component  of  the  outward  unit  normal  in  the  x. 
Xi 

tion,  is  applied  to  the  E terms.  This  yields 


3_ 

3x 


A 


g0 (r|r') 


_3 

3z 


Einc 


(r')ds' 


A 


gQ (*!*') 


3 2 ninc 
3x3z 


- S0(5|f')  3tEx"C(f')  Vr’>ds’ 

3A 


and 


gQ (*!*') 


_3 

3z 


EinC (r' )ds' 

y 


goC?|f)  4EynC(P:>  ny(*')ds'  • 

3A 


(4.18) 


direc- 


' )ds  ’ 


(4.19a) 


•)ds' 


(4.19b) 


Addition  of  the  first  two  terms  on  the  right  of  (4.19a)  and  (4.19b) 
yields 
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T(r)- 


JJ 


(flf)  i4Sx^“cs-)  * ^nC(f'))ds'  (4'20) 


g. 


Since  Einc  is  divergence-free  away  from  the  impressed  source.  Green's 
second  identity  (4.9)  gives 


T(?)»  - 


ijnc(?)  + 


3A 


irtf)  ^rg0c;if)]«' 


Thus  (4.17)  becomes 


_3 

3x 


g0(r|r')  ^x(r’)d£'  + ^ 


g0(r|r')  (r ' ) di ' 


3A 


3A 


= - 4 E^nC(?)  + 2 


tg0Ci|?'3 


3A 


(4.21) 


- 2 


f [^E“c(i'J  n y(r')  * \(f',lds' 


L3z  y 


3A 


3 

9y  JJ 


g (r|r' ) [T x (r' ) + T2(r’ ) ]dH’ 


3 
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g (r|r'J  [T-(r')  + Voids'  + 


Similarly  (4.18)  may  be  used  to  eliminate  the  derivatives  on  the 
integrals  involving  the  T Tj,  and  T4  terms  in  (4.22). 


In  conclusion  (4.13)  and  the  constraint  equation  (4.22)  enforced 
in  the  aperture,  plus  the  boundary  condition  that  the  tangential 
component  of  E vanish  on  the  perfectly  conducting  obstacle  (4.1c)  as 
well  as  on  the  aperture  boundary  is  sufficient  to  determine  the  electro- 
magnetic fields. 


CHAPTER  S 


NUMERICAL  FORMULATION  FOR  CAVITY  WITHOUT  OBSTACLE 


Numerical  solution  to  the  exterior-interior  coupling  problem 
of  Chapter  4 could  be  accomplished  by  discretizing  the  system  of 
integral  equations.  However  for  arbitrary  obstacle  geometries  the 
problem  requires  inordinate  matrix  sizes.  However  if  the  obstacle 
were  taken  to  be  a thin-wire  parallel  to  a coordinate  axis  only  one 
component  of  the  current  would  have  to  be  taken  into  account.  This 
would  reduce  the  matrix  size,  as  well  as  the  number  of  computations 
needed  to  solve  the  problem  discussed  in  Chapter  4.  Since  the 
numerical  difficulties  associated  with  the  integral  equation  obtained 
by  setting  the  tangential  § (4.1c)  to  zero  on  the  thin-wire  surface 
have  been  overcome  by  Seidel  [1977],  we  consider  the  exterior-interior 
coupling  problem  (Figure  1)  with  no  scattering  obstacle.  Equation 
set  (4.13)  for  the  aperture  fields  becomes 


{ [gpS’XX(r|r' ) + gc » xx ( - | - , ) ] EyC?,)ds, 


* 

gQ(r|r')  ^(r'JdA' 


* f - 

2 JJ  91  EjnC(r')ds' 


(5.1a) 


L_.l 


[gpS,/y (r|r')  + gp,yy(r|r')]  E (r')ds' 


= g0(r|r')  \jjx(f’)dJl' 


- 2 ff  g„ (?!*')  g|Ejnc(?')ds'  . 


(5.1b) 


The  constraint  equation  (4.22)  reduces  to 


2 J ir-EznC'f')  - Ez"C 

3A 

2 { g0(f|f,)  EynC ^ 


+ 3TExnC(?')  nx(f,)ldA’  • 


(5.2) 


In  the  above  equations  the  have  been  removed  from  the  E's  by 
multiplication  by  -iujy  (2.1c).  Similarly  the  auxiliary  unknowns  \p^  and 
ipy  have  been  replaced  by  and  \\>^  which  equal  -iwy  times  ^ and  ^ 
respectively.  The  electric  field  must  also  satisfy  the  boundary 
condition  that  its  tangential  component  vanish  along  the  perfectly 
conducting  aperture  edges. 
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The  method  of  moments  will  now  be  used  to  solve  (5.1)  subject 
to  (5.2)  plus  the  appropriate  boundary  conditions.  First  the  aperture 
fields  are  expanded  in  pulses. 


NX (NY-2)  1,  r e ith  x pulse 

E*(?)  = 1 exi  Pxitf)  PYiC?)  = { (5.3a) 


0,  otherwise 


NY (NX-2)  1,  re  ith  y pulse 

K&  = X ®yi  Pyi(?)  Pvi(f)  = f (S.3b) 


0,  otherwise 


where  NX  and  NY  denote  the  number  of  divisions  along  aperture  sides 
parallel  to  the  x and  y axes  respectively  (Figure  9) . In  general  the 
pulses  are  assumed  rectangular  in  shape.  Similarly  the  auxiliary  un- 
knowns are  expanded  in  boundary  pulses. 


V'3  = cxi  Psi(?)  V*3  = I cyi  Psi^ 


(5.4) 


where 


A,  r e ith  boundary  pulse 

Pgi(r)  = f } and  NT  = 2 (NX  + NY). 

0,  otherwise 


Next  the  system  of  integral  equations  is  converted  to  a matrix 
equation.  The  total  number  of  unknowns  in  (5.1)  and  (5.2)  is 
2NX*NY  + 2NX  + 2NY.  First  (5.1b)  is  enforced  at  nonzero  aperture 
pulse  midpoints.  Then  (5.1a)  is  enforced  at  nonzero  aperture  pulse 
midpoints.  This  yields  a total  of  2NX*NY  - 2NX  - 2NY  equations.  Next 
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(5.1b)  is  enforced  at  the  zero  side  pulse  midpoints  and  likewise  for 
(5.1a).  The  total  number  of  equations  is  now  2 NX*NY  for  2NX*NY  + 

2NX  + 2NY  unknowns.  It  is  sufficient  to  enforce  (5.2)  only  along  the 
aperture  boundary  [Mittra  et  al.,  1973].  Enforcing  (5.2)  at  the 
boundary  pulse  midpoints  yields  the  remaining  2NX  + 2NY  equations.  The 
resulting  matrix  equation  is  given  by 


GYY 

0 

GipX 

0 

~ 

e 

X 

FX“ 

0 

GXX 

0 

G^Y 

e 

y 

FY 

GYS 

0 

G4)XS 

0 

c 

X 

= 

FXS 

0 

GXS 

0 

G^YS 

c 

. y 

FYS 

_ 0 

0 

CX 

CY 

_FC  _ 

(5.5) 


The  partitioned  matrices  are  given  below. 


GYY(i.j)  = 


[g£S,yy(r  • |r')  + g£,yy(r  . |r')]ds' 
LBF  cxi1  SF  cxi1  1 


P . 


GXX(i.j)  = 


| [8SS,XXtfcyi|f')  * 4’XX(fcyilf,,lds' 


P . 
VI 


(5.6a) 


(5.6b) 


G^X(i, j ) 


G^Y(i.j) 


*o(*bi 


r')ds' 


(5.6c) 


(5. 6d) 


FXS(i)  = 2 g0(rxai|r’)  -gf  (i')ds' 


ff  9E 

FYS(i)  = - 2 jj  Cg0C?ysi|f-)  -£  (r')ds' 
A 


(5.7c) 


(5  • 7d) 


FC(i) 


- 4 E^nc(?b.)  ♦ 2 [g„(fM|r')  (?') 


5ov‘bi'*  ' 9n' 


9E  9E 

2 g0(fbil?')  CJ,)nyC?,5  if  C?')nxCr')]dI* 

9A 


(5.8e) 


The  major  obstacle  in  filling  the  needed  matrices  is  that  the 
cavity  dyads  in  (5.6a),  (5.6b),  (5.6e)  and  (5.6f)  given  by  Table  3 are 
slowly  converging  as  source  and  observation  point  become  close.  The 
summation  of  these  series  and  the  techniques  used  are  considered 
important  enough  to  merit  a seperate  chapter  (Chapter  6) . The  remaining 
matrix  elements  are  filled  first  by  analytical  integration  of  singular 
contributions,  then  by  Gaussian  quadrature  integration  of  the  remaining 
parts  [Abramowitz  and  Stegun  1970] . Evaluation  of  the  Hankel  function 
of  the  second  kind,  needed  for  g (r|r')»  is  done  by  means  of  a "poly- 
nomial approximation"  [Abramowitz  and  Stegun  1970] . To  evaluate  the 
forcing  functions  (5.7)  the  incident  fields  (4.14)  are  approximated  by 
the  centerpoint  values  on  each  aperture  pulse  (Figure  7).  For  normally 
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incident  excitation  this  approximation  becomes  exact.  The  only  other 
obstacle  to  matrix  fill  is  the  bookkeeping  necessary  to  insure  that  all 
elements  are  computed  with  the  least  amount  of  computations  as  possible. 

The  matrix  of  the  linear  system  (5.5)  can  be  numerically 
inverted  in  a manner  that  makes  advantage  of  the  large  blocks  of  zeros. 
Once  this  matrix  is  inverted  it  may  be  used  to  solve  the  aperture  field 
responses  due  to  different  polarizations  and  incidence  angles  of  the 
incident  plane  wave  field.  The  method  of  partitioning  is  used  to  invert 


the  matrix. 


It  is  straightforward  to  show  [Faddeev  and  Faddeeva  1963]  that  if 


C J D 


(5.9) 


where  A and  D square  matrices  of  order  p and  q respectively,  then 


K ! L 


(5.10) 


where  K and  N are  square  matrices  of  order  p and  q respectively.  The 
partitioned  matrices  of  S * are  given  by 


N = (D  - CA-1B) -1 


M = - NCA 


L = - A-1  BN 


K = A"1  - A-1BM 


(5.11) 


provided  that  each  inverse  exists. 


ik 


This  scheme  can  be  applied  to  the  matrix  of  (5.5)  where 


A = 

GYY  ! 0 ' 

4. 

. 0 ] GXX_ 

B = 

"g<|jX  ! 0 ' 

. 0 ! Gt|)Y_ 

GYS 

o ' 

*G1|)XS  ! 0 1 

C = 

o 

GXS 

D = 

0 

1 G^YS 

1 

O 

— 

o 

J 

cx 

fcY 

GYY  and  GXX  are  square  matrices  of  order  NX(NY-2)  and  NY(NX-2).  For 
large  NX  and  NY  these  matrices  will  dominate  considerations  of  computer 
storage.  Matrices  B,  C,  and  D are  [NX(NY-2)  + NY(NX-2)]  x 4(NY  + NY), 

4 (NX  + NY)  x [NX (NY- 2)  + NY (NX-2)],  and  square  4(NX  + NY)  matrices 
respectively.  Once  K,  L,  M and  N are  found  the  fields  can  be  obtained 
by  matrix  multiplications. 


CHAPTER  6 


EVALUATION  OF  CAVITY  POTENTIAL  TERMS 


The  final  obstacle  to  solution  of  the  exterior-interior  coupling 
problem  of  Chapter  5 is  the  matrix  fill  of  elements  involving  the  cavity 
vector  potential  gp  in  (5.6).  Table  3 gives  this  dyad  represented  as  a 
threefold  infinite  eigenfunction  expansion.  For  source  and  observation 
points  distant,  it  will  be  seen  that  this  threefold  infinite  eigen- 
function expansion  reduces  to  a rapidly  exponentially  convergent  double 
series.  However  as  source  and  observation  points  become  close,  the 
series  becomes  slowly  convergent.  This  chapter  develops  an  efficient 
means  of  evaluation  of  the  cavity  potential  dyad  for  source  and 

r 

observation  point  close  under  the  restriction  that  ka  (Figure  1)  is  less 


than  it. 


Consideration  is  now  given  to  the  xx  component  of  the  cavity 


C XX  c w 

dyad  g ’ . Similiar  techniques  hold  for  g ’ . From  Table  3 with 

r 


z = z ' = c (Figure  1) 


gp’XX(?|?') 


°°  £. e e 

-r—  1 ^ n sin(k  x)sin(k  x')cos(k  y)cos(k  y') 

abc  1=0  K»  -k2  x r Yy 

, x.mn 
m=l 

n=0  (6.1a) 


= S ( | x-x' | , | y-y ' | ) + s (|x-x'|,  y+y')  - S (x+x',  y+y') 


- S (x+x',  | y-y' I) 


(6.1b) 
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where 


stx’Y>  - IffiT  j„  “>»»/)  ' 

, KS.mn  “ k 
m=l 

n=0 


(6.1c) 


Since  the  terms  of  (6. lb)  depend  only  on  sums  and  differences 
of  the  respective  source  and  observation  points'  components,  the  number 
of  calculations  in  evaluating  (5.6b)  and  (5.6f)  is  less  than  if  (6.1a) 
were  applied  directly. 

For  source  and  observation  points  distant  the  threefold  infinite 
series  of  (6.1c)  may  be  reduced  to  a double  sum  that  yields  itself  to 
rapid  evaluation  by  digital  computer.  Equation  (6.1c)  is  closed  in  Y 
[Gradshteyn  and  Ryzhik  1965]  which  yields 


. 00  cosh[A2  + k2  - k2  (b-Y)] 

SCX.Y)  = J-  l 4 —5 7 — cosCk  X) 

Aac  1=0  A2  + X2  - k2  sinh  [A2  + k2  - k2  b]  x 

.XV  XV 

m=l  1 J 

(6.2a) 


with  ka  (Figure  1)  less  than  rr.  The  series  in  (6.2a)  is  exponentially 
converging  and  readily  evaluated  provided  that  Y is  not  too  small  with 
respect  to  the  aperture  dimensions  and  not  too  close  to  2b.  Likewise 
if  (6.1c)  is  closed  in  X one  obtains 


S(X,Y)  - ^ l0  T(X)  cosoyn 

n=0 


1 y £&£ncosCkyY^ 

4abc  H=0  k2  + k2  - k2 
n=o  y 2 


(6.2b) 
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where 


if  k2  + k2  - k2  > 0 

y z 


if  k2  + k2  - k2<  0 

y z 


Similar  techniques  may  be  used  to  close  the  second  term  of  (6.2b).  The 
result  is  a single  infinite  series  with  known  asymptotic  form 
[Abramowitz  and  Stegun  1970,  p.  1005].  Thus  this  series  may  readily  be 
evaluated.  First  termwise  subtraction  of  the  corresponding  asymptotic 
term  from  the  original  series  term  yields  a rapidly  converging  new 
series.  Next  the  otherwise  slowly  converging  asymptotic  series  is 
added  to  the  modified  series  by  use  of  its  known  analytical  form.  The 
first  term  of  (6.2b)  is  exponentially  convergent  provided  X is  not  too 
small  or  too  close  to  2a.  In  the  aperture  geometries  considered 
(Figure  1),  the  only  case  where  neither  (6.2a)  or  (6.2b)  yields  itself 
readily  to  numerical  evaluation  is  for  both  X and  Y small  compared  to 
the  cavity  dimensions.  Otherwise,  Seidel  [1977]  has  obtained  an 
efficient  method  for  numerical  evaluation  of  these  sums. 

The  remainder  of  this  chapter  is  concerned  with  evaluation  of 
(6.1c)  near  the  source,  that  is  for  X = |x-x'|  and  Y = | y-y ' | small. 
Equation  (6.1c)  is  closed  in  z,  making  use  of  Spiegel  [1964,  p.  189] 
and  Gradshteyn  and  Ryzhik  [1965,  p.  40].  This  yields 


—7T7T 
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CO 

S Cx» Y)  = 2ib  I en  T(k2  + k2  - k2)  cos(kxX)  cos(kyY)  (6.3) 


where 


T<kx  + 


, if  k2  + k2  - k2  > 0 

a y 


if  k2  + k2  - k2  < 0 . 
x y 


This  series  in  (6.3)  has  an  asymptotic  form 


WX’Y> 


CO 

= — y 

2ab  L 


n=0 


1 A 


k2  - 

y 


cos(k  X)  cos(k  Y) 
.2  x y 


(6.4) 


The  difference  series  whose  terms  consists  of  the  original  series  terms 
(6.3)  minus  the  corresponding  asymptotic  series  terms  is  readily  evalu- 
ated by  digital  computer.  Thus  (6.3)  may  be  evaluated  if  a method  of 
evaluating  the  otherwise  slowly  converging  asymptotic  series  may  be 
found. 

To  evaluate  SASy(X,Y)  (6.4)  consider  the  single  series 

So°°  = X cos(k  Y)  (6.5) 

0 n=0  A2  + u2  y 

where  Y is  greater  than  zero  and  u2  equals  k2  - k2  which  is  greater  than 
zero  under  the  restriction  that  ka  is  less  than  it.  An  alternate 
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expression  for  Sq(Y)  will  be  obtained  from  Poisson's  formula  [Papoulis 


1962] 


l f(nT)  = y l F(^r> 

n=-<r 


(6.6a) 


where  f(t)  is  continuous  and  its  Fourier  transform  F(w)  is  given  by 


F(w) 


■r 


f(t)e'iwtdt  . 


(6.6b) 


Equation  (6.5)  is  rewritten  as 


ttY, 


S (Y)  = Y l C°S(n  b} 


n=-°°  / „_v  2 2 
(MX)  +V 


(6.7) 


where  v equals  Yu.  In  (6.6a)  f(t)  is  identified  as  ' — . Its 


/ 1 2 + v2 


Fourier  transform  is  given  by 


F(u)  = 


f00  -1U)T 


(6.8a) 


f cos[(l  ♦ m)tl  dt  + 

A2  ♦ v2 


cos  r (ai-l)tl 


-'T2  + v2 


(6.8b) 


Kq  ( 1 1 + <u|v)  + Ko(|u  - l|v) 


(6.8c) 


j 


80 


Equation  (6.8b)  follows  from  (6.8a)  plus  symmetry  considerations.  The 
integrals  of  (6.8b)  result  in  the  modified  Bessel  functions  of  (6.8c) 
[Morse  and  Feshbach  1953,  p.  1323].  Application  of  (6.6a)  to  (6.7) 
yields 


VY>  ■ ? I EK0C|y  ♦ 


cj  - k2)  ♦ Ko(|2nb  - Y | 


c2  - k2)] 
x J 


(6.9a) 


^ Kq(^2  - k2 Y)  + ^ l [Kn([2nb  + Y]A2  - k2) 


IT  ^,0 
n=l 


♦ KQ([2nb  - Y] A2  - k2)] 


Application  of  (6.9b)  to  (6.4)  yields 


(6.9b) 


w ■ ....  - 

S4oY(XA)  - ±=  l K (A2  - k2  Y)  cos(kX) 


air  o'  x 

m=l 


+ h l l CKoa2nb  + Y]A2  - k2) 


air  u.  o 

m=l  n=l 


+ K ([2nb  - Y] A2  - k2)]  cos(k  X)  . 

OX  X 


(6.10) 


The  second  term  on  the  right  of  (6.10)  is  very  rapidly  convergent  due 
to  the  exponential  behavior  of  Kq's  asymptotic  expansion  [Arfken  1971, 
p.  517].  The  first  term  on  the  right  of  (6.10)  is  exponentially 
convergent  for  Y > 0.  However,  as  Y becomes  smaller  the  series  converges 


Note  that  the  dominant  singularity  in  (6.11)  is  the  same  as  that  of  the 
half  space  dyadic  gp  of  (3.9a).  Unfortunately  the  infinite  series 
(6.11a)  is  slowly  convergent.  However  if  one  approximates  each  term  of 
this  series  by  a finite  Taylor  series  in  X and  Y about  the  origin,  the 
summation  over  n may  be  obtained  by  analytical  methods.  For  example, 
the  first  term  obtained  in  the  Taylor  series  is  given  by 

l -~a~an)-  = - *"[2  sin(ka)]  (6.12) 

n=l 

[Abramowitz  and  Stegun  1970,  p.  1005].  Similar  formulas  may  be  used  to 
evaluate  the  higher  order  terms  in  the  Taylor  series  approximation. 

Although  the  above  discussion  has  been  restricted  to  the 

c xx  c yy 

evaluation  of  gp’  only  minor  changes  are  needed  to  evaluate  gp  /. 

Also  the  terms  that  must  be  evaluated  (5.6)  are  integrals  of  the  xx  and 

y?  components  of  g£.  The  sums  of  (6.2)  may  be  evaluated  by  term  by  term 
F 
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integration  [Titchmarsh  1939,  pp.  36-45].  The  exponentially  convergent 
terras  (6.10)  are  also  evaluated  termwise.  The  only  remaining  question 
is  the  accuracy  of  the  previously  discussed  Taylor  series  approximation 
(6.11).  A sixth  order  Taylor  series  was  used  for  numerical  testing  and 
yielded  no  relative  error  worse  than  1*10  ^ in  these  terms  when 
max(X,Y)/a  < .05.  This  procedure  allows  rapid  evaluation  of  the  series 
on  the  right  of  (6.11),  provided  that  both  X and  Y are  small  with  respect 
to  the  cavity  dimensions  in  the  aperture  plane.  Otherwise  this  series 
would  not  readily  yield  itself  to  numerical  evaluation  on  a digital 
computer . 


CHAPTER  7 


CONCLUSION 

In  this  paper,  the  task  was  undertaken  to  formulate  the  problem 
of  an  aperture  in  a perfectly  conducting  infinite  plane  backed  by  a 
rectangular  cavity  with  a wire  obstacle.  The  work  of  Seidel  [1977, 
1978]  was  extended  to  include  apertures  large  enough  so  that  the  dipole 
moment  approximation  is  no  longer  necessary.  Although  the  theory  pre- 
sented is  valid  for  more  general  shapes  of  obstacles,  the  particular 
problem  of  interest  is  a thin  wire  parallel  to  a coordinate  axis.  This 
will  reduce  matrix  sizes  and  computation  time. 

One  essential  aspect  of  this  formulation  has  been  the  dyadic 
Green's  functions.  Their  completeness  properties  and  singularities 
received  special  attention.  Another  important  result  was  the  summation 
technique  of  Chapter  6 under  the  restriction  that  ka  (Figure  1)  is  less 
than  tt.  This  technique  allowed  summation  of  magnetic  source  vector 
potential  dyads  with  source  ^nd  observation  poin:  , close.  Without  this 
technique  the  dyadic  series  would  become  very  slowly  convergent  as 
source  and  observation  points  become  close.  This  result  made  consider- 
ation of  large  apertures  and  the  removal  of  dipole  approximations 
possible. 

Future  Work 

At  this  time  numerical  results  should  be  possible  for  the 
problem  without  the  wire  obstacle  (Chapter  5).  The  computer  program 
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should  be  checked  by  removal  of  cavity  and  comparison  with  the  litera- 
ture, especially  the  case  of  small  apertures.  The  problem  of  large 
matrix  sizes  and  large  computer  run  times  could  be  reduced  by  symmetry 
considerations . 

Addition  of  an  obstacle  of  general  shape  would  greatly  increase 
matrix  sizes.  However  if  a thin  wire  parallel  to  two  cavity  sides  were 
used,  the  problem  should  be  feasible.  The  major  numerical  problems 
involved  in  addition  to  this  thin-wire  have  been  overcome  by 
Seidel  [1977,  1978].  The  major  remaining  work  is  computation  of 
coupling  terms  between  wire  and  aperture  fields.  In  addition  to 
symmetry  considerations,  the  use  of  a slot  with  the  appropriate  edge 
behavior  of  the  fields  built  into  the  formulation  (Chapter  3)  would 
enable  one  to  consider  a one  dimensional  array  of  unknown  aperture 
field  elements.  This  procedure  may  somewhat  complicate  numerical 
considerations,  but  it  would  greatly  reduce  matrix  sizes. 

Finally  the  transient  problem  needs  to  be  considered.  The  time 
domain  solution  can  be  obtained  by  inverse  Fourier  transformation.  It 
is  noted  that  many  frequency  values  are  needed.  Some  type  of  symmetry 
considerations  or  the  previously  discussed  slot  considerations  should  be 
considered  to  reduce  computation  time. 


APPENDIX  A 


SOME  SYMMETRY  PROPERTIES  OF  DYADIC  GREEN'S  FUNCTIONS 

The  symmetry  properties  for  and  gh  given  by  (2.9)  and  (2.14) 
are  now  verified.  Consider  Equation  (2.6)  with 

P = a • Se(r|r')  , 

where  a is  an  arbitrary  constant  vector  and 
Q = 8e(r|r") 

yields 

ffr  _ r _ 

jjJ  * Ge(r|?')  * V x V x Ge(r|r")  - (V  x V x a • 5e(r|r')) 
v 

• 5e(r|r")]dv  = 0 , 

since  the  surface  integral  vanishes  by  (2.11b)  and  the  radiation  condi- 
tion. Noting 

a • Sg(r |r")  = §e(f|f")  • a , 

and  using  (2.7)  yields 

a • G (r"|r')  = a • 8 (r ' |r")  . 

“ G 
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one  obtains 
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[(V  x V x Ggfrlfp  - k2§e(r|r1))  • a • V x gh(r|r2)  • b 

J J 

v 


- V x 5g(r|r1)  • a • (V  x 7 x Ih(r|r2)  - k2gh(r|?2))  • b]dv 


[§e(r|r1)  • a • V x gh(r|r2)  • b 


- V x ^(rlr^  • a • gh(r|r2)  • b]dv 


||  [(V  x Se(r|r1)  • a)  x (V  x gh(r|r2)  • b)]  • n ds 


The  surface  term  vanishes  for  r^  and  r9  interior  to  v by  the 
radiation  condition  plus  the  boundary  condition  that  n x gg(r|r2)  = 0 
on  a perfectly  conducting  surface.  By  (2.7c)  and  (2.12c)  one  obtains 

£ • V1  x Ih(?1|r2)  • 5 - b • V2  x Sg(r2 |rp  • £ 

= - k2  HI  [§e(r|r1)  • a • V x gh(r|?2)  • b 
v 

- V x 5e(r|r1)  • £ • gh(r|r2)  • b]dv  . 

Once  the  volume  integral  is  shown  to  be  zero  we  will  have  our  final 


result 
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[V2  x = V1  x • 

By  the  divergence  theorem  the  volume  integral  becomes 

" k2  II  [(gh(r|?2)  • b)  x Se(r|r1)  • a]  * n ds 
S 

which  is  zero  for  r 2 and  r^  interior  to  v since 
n x 8e(r|r1)  = 0 

on  a perfectly  conducting  surface  and  the  radiation  condition  at 
infinity  holds. 
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